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ABSTRACT 



Based on the theory of continued fractions, a technique 
is developed for the reduction of high order multivariable 
systems. The mathematical basis for which these techniques 
work is elucidated, and its superiority of the mixed form 
over any other form of continued fractions is established. 

The general solution to linear regulator problem is developed 
and the properties which this solution exhibit are elucidated. 
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I . INTRODUCTION 



The capability of obtaining simplified mathematical 
models for use in the analysis of high order dynamic systems 
has traditionally relied on the experience and ingenuity of 
the analyst. Usually, these efforts have been achieved 
using both frequency and time dcmain techniques. 

In dealing with the problem stated above, and through 
the use of theory of continued fraction, this research has 
found a series of properties applicable to reduction of 
multivariable systems, lower order observers and derivations 
of lower order systems for the linear regulator problem. 

The relationships developed, show applications in areas 
where the high order systems are impractical or undesirable 
to use due to their complexity or difficulty of implementarion. 

In Chapter II reduced order models of multivariable 
systems for the first, second and mixed Cauer forms are 
developed. Techniques for approximating a high order linear 
rime-invariant system with various inputs and various outputs 
by a reduced order model, have been suggested by Chen Cl], 
Meier, L and Luenberger [2], L. S. Shieh and Y. J Wei [3], 

M. R. Calfe and M. Healey [4], L. S. Shieh, J. M. Navarro 
and R. Yates [5], D. A. Wilson [6], L. S. Shieh and F. F. 
Gaudiano [7]. 
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Most of these methods for reducing high order linear 
systems are based on the following principles: 

1. The low performance terms can be discarded and the 
high performance terms should be retained. 

2 . Linear transformation to obtain matrix diagonaliza- 
tion where certain diagonal elements can be dis- 
carded . 

3. The sum of squares of the errors between the 
responses of the real system and those of zhe 
approximate model at the sampling instant is 
minimized in order to obtain the parameters of the 
approximate model. 

In a recent paper, Chen [l] proposed a reduction of 
multivariable control systems by means of matrix second 
Cauer form of continued fractions. Through the method, a 
simplified model is obtained by keeping the first several 
significant matrix quotients and discarding the others. 
However, the technique (due to the nature of the Cauer 
second form) , provides satisfactory results in the steady 
state region only. Furthermore, M. R. Calfe and M. Healey 
[4], have shown that the method does not guarantee the 
reduced model to be stable. 

In Chapter III derivation of lower order system for the 
linear regulator problem via Cauer form is obtained and 
also a near optimal solution for the original system can be 
found through a reduced system. 
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II. MULTIVARIABLE SYSTEMS REDUCTION VIA THE CAUER FORM 



A. THREE MATRIX CAUER FORMS 

L.S. Shieh and F.F. Gaudiano [7] have shown that in terms 
of multivariable systems the quotients in the three Cauer 
forms are replaced by matrix quotients and the division in 
the continued fraction is replaced by matrix inversion. The 
first matrix Cauer form is 



T ( s ) = [H’s + [H^ + [H' 3 s + [H^ + C . . . I" 1 ]" 1 ] -1 ] -1 ] 1 (1) 

the second matrix Cauer form is 

T(s) = [H 1 + [H i + [H, + [H i + [ . . .]" 1 ]" 1 ]" 1 ] -1 ]" 1 (2) 
X 2s 3 4s 

and the mixed matrix Cauer form is 



T(s) = [K.. + KLs + [K_— + K’ + [K_ + K's + K„- 
1 1 2s 2 3 3 4s 

+ K* 4 + C . . .]“ 1 ]" 1 ]" 1 ]“ 1 ]" 1 (3) 

where H^, H^, , and K'^ are constant m x m matrix quotients 

obtained respectively from the matrix Routh's array and the 
generalized Routh’s algorithm shown in equations (4a) and 
(4b) . 



1 1 



H n 



H, 



H, 



A 11 A 21 



A A ’ 1 
21 31 





11 12 



21 22 



A 



A 31 A 41 






31 32 



A, A, -i 

In l,n+l 



A 



2n 







HU = 



41 42 



'3 ,n-l 



4 ,n-2 



^>*3 = 



A. , A “ 1 
1 ,n+l 2n 



A 0 A 0 , 
2,n 3,n-l 



A A " 1 

3 , n-1 4 ,n-2 



H , = A , 1 A 1 . ( 
n-1 n-1 ,1 n,l \ 



A n-1 , 1 A n-1 , 2 A n-1 , 3 



H = A .A 1 .. . 
n n,l n+1,1 



n , 1 n , 2 



n+1,1 



> 




H* , =A , 0 A _1 - 
n-1 n-1, 3 n,2 



H=A 0 A 1 , . 
n n , 2 n+1,1 



(4a) 



K„ 



K, 



K, 



A 11 A 21 



A 21 A 31 



i< 



A 11 A 12 



At A- ,, 
In l,n+l 



A 21 A 22 



< 



A 31 A 41 






A 31 A 3 2 



A 41 A 42 



2n 



3 ,n-l 



4 ,n-2 



> 

> 



3 * 



KL = 



K' = 



A. .A ' 1 
l,n 1 2n 



A 0 A . 1 
2 ,n 3 ,n-l 



A A -1 

3 ,n-l 4 ,n-2 
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K , =A . , 1A" 1 
n-1 n-1 5 n 



A , , A -jo 
n-1,1 n-1,2 



-1,1 < 



A , A 0 
n , 1 n , 2 



K =A , 1A , 
n n 5 n+1,1 



< 

A 



n+1,1 



X n-l =A n-l,3 A n-2 

> K n =A n,2 A n+l,l 



(4b) 



where the elements of the first and second rows of equations 
(4a) and (4b) are the matrix coefficients of the system 
given in equation (5). 



T ( s ) = [A 2 ,n S n_1 + A 2 n-lS n " 2 +...A 23 S 2 + A^S + A^jx 



(5) 



[A l5 n+lS n + A x ,nS n_1 +...A 13 S 2 + A l2 S + A 11 ]“ 1 



The elements of the third, fourth and subsequent rows in (4a) 
and (4b) are evaluated respectively for the three Cauer forms 
by the formulation shown in (6a) through (6c). 

Hp = Ap,n+2-pAp^ +1,1 p = l,2,3,...n 

T 

Aj,i = Aj-2, i+1 - Hj-2 Aj-l,i j = 3,4, ...2n 

det Ap+1,1 ^0 i = n+1, n, n-l,...l 

(6a) 



13 



Hp = Ap , lAp" +1,1 



p = l , 2 ,3 , . . . 2n 



Aj,i = Aj-2, i+1 - Hj-2Aj-l, 1+1 j=3,4,...2n 
det Ap+1,1 i 0 i=l,2,3,... 



(6b) 




p = l , 2 , 3 , . . . n 



j=3,4,.. .n+1 



A j , i = Aj-2, i+1 - Kj-2Aj-l, i+l-K* j-2 Aj=l,i 1=1, 2, 3, . 



(6c) 



It is important to note that since the Cauer first and second 
forms are special cases of the Cauer mixed form, their 
formulation in (6a) and (6b) can be derived directly from 
(6c), by letting all kp's or all K£> ' s equal to zero respec- 
tively. 

B. STATE SPACE FORMULATION FOR THREE CAUER FORMS 

The Cauer Mixed Form - Consider a typical feedback system 
with a minor feed forward loop as shown in Figure 1. The 
overall transfer function is given by: 



Equation (7) can be rewritten as a mixed matrix of continued 
fractions 



T ( s ) = = [G + F] [I + (G + F)H] 1 . 



(7) 



T ( s ) = [H + [F + G] 1 ]“ 1 



( 8 ) 
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where 



H = K x + K 1 S \ 

F = K 0 — + Kb ( (9) 

2s 2 l 

I = Identity Matrix' 



U( s ) 




Y(s) 



Figure 1. Feedback and Feedforward controls. 

If the subsystem G is expanded again, equation (3) is obtained. 
This equation can be represented by the block diagram shown 
in Figure 2. Where a 2 inputs-2 outputs nth order system 
is shown. Again, it is important to note that if all Ki's 
go to zero in Figure 2 the block diagram representation of 
Cauer matrix first form nth/ 2 order system as shown in 
Figure 3 will automatically be obtained. In a similar 
fashion if the k’i's go to zero, the block diagram represen- 
tation of Cauer matrix second form nth/ 2 order system as 
shown in Figure 4 is obtained. 
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Figure 2. General Matrix representation of a nth order system by 
Cauer mixed form with two inputs two outputs, (n even) 
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Figure 4. Matrix representation of a nth/2 order system by Cauer 
second form with two inputs two outputs, (n even) 
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Going back to Figure 2 and allowing to be the state 
variable vector (same order as the matrix K or K' ) , the 
equations in the time domain can be written as follows: 



K i K 2"l + K i% g 2 + K i K 6^3 + ••• K i K n*n/2 



-(I + + K 1 K’ 2 )e 1 -(K' 1 K 4 + ) e 2 - (K^Kg+K^g ) e 3 



• (K iV¥n^n/2 ' K l K 2 e l “ K l K 4 e 2 - * * - " K l K n e n/ 2 + 



IU 



( 10 ) 



K' Kle, + ( K 1 , +Kt ) Kf e 0 + (K 1 , +K!L ) KLe - + • . . + (KI +K* )K* e /0 
121 1342 1363 13nn/2 



-(K 1 K’ 2 + K , 1 K 2 )e 1 -[(I + K! l +K' 3 )K 4 +(K 1 +K 3 )K’ 4 ]4 2 - 



[ (KL +Ki )K-+(K, +K q )K' k ]e 0 - . . . - [ (K*, +K~ )K +(K 1 +K Q )K' ]e /0 
13b 13b3 1 3 n 13nn/2 



-K,K 0 e,-(K,+K Q )K.e 0 -(K, + K 0 )K C -. . .-(K,+K_)K e ,_ + IU 
121 1342 136 13nn/2 



Ki K ^l +( H + 4 ) ¥2 +( ¥^3^ 5 ) K^3 + - + ( ¥K3 +K 5 ) ¥n/2 



- ( K 1 K , 2 +K' 1 k 2 )4 1 - [ ( Ktj+K’g ) K 4 + ( K 1 + K 3 ) ] 4 2 - [ 1+ ( Kt^+K^ + KV ) K g + 

( V K 3 + K 5 )K 6 ]4 3---- C(K ¥ K 3 + K 5 )K n +(K l + K 3 + K 5 )K, n ] ¥/2 



( 11 ) 



~ K l K 2 e l _(K l + K 3 ^ K 4 e 2 < ‘ i< l +K 3 + K 5 )K 6 e 3 (K 1 + K 3 + K 5^ K n e n/2 + IU 



( 12 ) 
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K i K, 2 S l +(K i +K 3 )K, 4 S 2 +<K i +K, 3 +K 5 >K, 6 S 3 + •” 



+ (K' 1+ K' + K' 5+ ... + K^ 1 )K^ n/2 



-CK' K 2 +K 1 K^)e 1 - [(K| L +K , 3 )K 4 +(K 1 +:< 3 )K' 4 ]e 2 



[(K'.+K' +K' t -)K,t(K.tK,tK c )KV]e,-. . 
1356 135o3 



CI+(K T ,+K' +K T ,+ . . .+K T , )K + (K 1 +K,+K I -+. . . +K , )K' ]e /0 
13 5 n-1 n 13 5 n-1 n n/2 



-K 1 K 2 ei-(Ki + K 3 )K 1(e 2- ( V K 3 +K 5 ,K 6 e 3- 



(K 1 +K 3 + V ••• * K n-l )K n e n/2 + IU 



and 



(13) 



Y = K 2 e 1+ K' 2 4 1+ K 4V K' 4 e 2+ ... +Vn/2 * K ’ n e n/2 



(14) 



where U and Y are the input and output variables of the 
system respectively and I is the Identity Matrix. From 
Equations (10), (11), (12), (13), and (14) the following 
corresponding matrix formulation can be obtained: 

(15) 

(16) 



[A 1 ]E = -CA 2 ]E - [A 3 ]E + [ 3 ] U , 

and 

Y = [C 1 ] T E + [C 2 ] T E. 
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where 



[A 1 ] = 



KJK^ 

K 1 K 2 (k [ +K 3 )K 4 



K i K 6 



' K 1 K A 



(K’+K' )K' (K’ +K' )K' 

13 6 1 3 n 



K 1 K 2 (K i + k ’ 3 )K 4 (K i +K 3 + K 5 )K 6 



(K’ +K’ +K’ )K’ 
1 3 b n 



K1 K’ ( K', +K1 )K! (K* +K' +K1 )K’ . . (K’ +K' + . . . +K’ . )K' 

12 1341356 13 n-1 n 



[E ] = 



'n/2 



[a 2 ] = 



I + K»_K 2 +K 1 K , 2 


K 1 K 4 + K i K 4 • * ‘ 




K 1 K 2 +K 1 K 2 


1 + (K 1 + K 3 )K 4 + 


(KJ_ + K’ 3 )K 


K 1 K’ 2 + K! l K 2 


(K 1 + K 3 )K 4 + (K i 


+ K’ 3 ) K 4 . . . 


_ K l K 2 ;K i K 2 


(K 1 + K 3 )K 4 + (1 ^i 


+ K’ 3 )K 4 . . . 


. K, K’ +K’ K. . . . 
lb 16 







(K 1 +K 3 )K 6 +(K 1 +K 3 )K 6-*- 

I+(K 1 +K 3 +K 5 )K 6 +(K 1 +K 3 +K 5 )K 6--- 



(17) 



(18) 



(K 1 +K 3 +K 5 )K, 6 + (K 1 + K 3 + K 5 )K 6-*- 
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. . . K.. K f +K'K 
In In 

. . . (K;+KL)K + (K n +K-)K' 

1 3 n 1 3 n 

. . . (K* +K' +K')K + ( K, +K 0 +K c ) K* 
1 3 S n 1 3 d n 



. .I+CK’+K’ +K’+. . .+K' , )K +(K n +K 0 + K c +. . .+K ,)K r 

13 5 n-1 n 13 5 n-1 n 



(19) 



HE] 



( 20 ) 



'n/2 



K 1 K 2 K 1 K 4 



K 1 K 6 



K 1 K 2 (K 1 +K 3 )K 4 (K 1 +K 3 )K 6 
[A 3 ]= K x K 2 (K 1 +K 3 )K t+ (K 1 +K 3 + K 5 )Kg 



K !K 2 (V K 3 )K 4 (K 1 +K 3 +K 5 )K 6 



K n K 

1 n 

. (K.. +K„ )K 
1 3 n 

. (K n +K-+K, )K 
1 3 5 n 



( K, +K 0 + K C + • • • +K , )K 
135 n-1 n 



( 21 ) 
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e 



1 



e 



2 



[E ] = 



e 



3 



e n/2 



I 



I 



[ B ] = 



I 



I 



[C 1 ] T = CK 2 K 4 Kg ... K n ] , 



[C 2 ] T = [K * 2 K > 4 Kg ... K^] , 

















— ~ 




*1 




x m+l 




X 2m+1 




X (n/2-l)m+l 


e l 


*2 


e 2 


X m+ 2 


g 3 = 


x 2m+ 2 


®n/2 = 


X (n/2-l)ra+2 




X m 




X 2m 




X 3m 




X (n/2)m 



( 22 ) 



(23) 

(24) 

(25) 

(26) 
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*1 




x m+l 




X 2m+1 




x (n/2-l)m+l 


X 2 




X m+2 




X 2m+2 




X (n/2-l)m+2 


# 

# 

X 

_ m 


*2 = 


_ X 2m_ 


4 3 = 


♦ 

• 

_ x 3m 


e n/2 = 


_ x (n/2)m _ 



(27) 





X 1 




x m+l 




X 2m+1 




X (n/2-l)m+l 


e l = 


X 2 

• 

• 


e 2 = 


X m+ 2 


e 3 = 


X 2m+2 


e n/2 = 


x (n/2-l)m+2 




X 

__ 




_ X 2m_ 




CO 

B 




_ x (n/2)m _ 



(28) 

where m is the number of inputs and outputs. It is important 
to emphasize the properties that the Cauer matrix mixed form 
exhibits. If again the Ki's go to zero in the state formula- 
tion described by equation (15) and (16), they will reduce to: 

[A ]E = -CA 21 ]E + [B]U , (29) 

and 

Y = [C 21 ] T E. (30) 

and letting K'i's = H’i's equation (29) and (30) define the 
state space formulation of the Cauer matrix first form for 
an nth/ 2 order system, where n is even, Equations (29) and 
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(30) are second order differential equations which can be 
simplified to first order differential equations by 
assigning a new state variable [Z] = [E] . Thus equations 
(29) and (30) can be rewritten as: 



CA 11 ]Z ■ 


Ca 21 ]z + 


[3 ]U , 






and 












Y = 


[C 


] T Z. 
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where [A.^] 


= lim [A^] 


and K ’ i ' s 


= H'i’s 








Ki’s -*■ 0 




_ 




then : 


H i H 2 


T H, 4 


H 1 H 6 


. . . H' H’ 
1 n 






H 1 H 2 


(H' 1+ H')H' 4 


(Hi +H’ )H’ 
lob 


. . . ( h'^+h' 3 


> H ’n 


II 

1 — 1 
i — 1 

rH 

< 

i — 1 


H 1 H 2 


(H' 1+H ' 3 )h' 4 


(H i +H 3 +H 5 


)H 6-*- (H i +H 3 


+H' 

5 n 




H 1 H 2 


(H'+H', )H' 


(H i +H 3 +H 5 


)H’ . . . (HI +H’ 
6 13 


+H 5 + -* 



I' , )# 
n-1 n 



Z = 



'n/2 



( 31 ) 



(32) 



(33) 



( 34 ) 
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where e^, e 3 ... e n / 2 are S iven in equation (27) 



[A 2 i] = lira CA 2 ] = [I], 
Ki ' s-»-0 



( 35 ) 



where I is the Identity matrix, 



'n/2 



(36) 



where e^, e 2 , e^ ... e n / 2 are gi ven i- n equation (28) 

[C 21 ] T = lim[C 2 ] T and K'i’s = H’i's [C 21 ] T = [H' 2 H' 4 H'g-.H^] 



Ki's+O 



(37) 



This formulation given by equations (30) through (37) could 
have been obtained directly by inspection from Figure 3 . 

Similarly, if all K'i's go to zero in equations (15) and 
(16) then 

E = -CA 32 ]E + [ B ]U , (38) 

and 

Y = [C 12 ] T E, ' (39) 
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where 



[A 32 ] = liraCAg] and Ki's = Hi's 
K'i's+O 



[a 32 ]= 



H 1 H 2 H 1 H 4 



H 1 H 6 



. . .H.H 
1 n 



H, H 0 (Ht+H^JH., (H,+H q )H c . ..(H,+H,)H 
1 2 1 2 4 loo ion 

H 1 H 2 (H 1 +H 3 )H 4 (H l +H 3 tH 5 )H 6"' (H l +H 3 +H 5 



)H 



n 



H, H 0 (H +H„ )H. (H,+H,+H,)H C ...(H 1 +H,+H,+...K , )H 

12 1 O r 1 O b O lob n-l n 



(40) 



and 



[C l2 ] i = limEC^] 1 and Ki's = Hi's 
K ' i ' s-*0 



Cc i2 ] = [h 2 h 4 H 6 ... H n 3. 



(41) 



which define the state space formulation for the Cauer matrix 
second form. Again, this preceding matrix formulation could 
have been derived directly by inspection from the block dia- 
gram shown in Figure 4. Equations (40) and (41) have the 
same form of the state space formulation given by Chen [1]. 

At this point, the following observations may be made re- 
garding the state space formulation just developed: 

1. The .elements in the state matrix [Ai’s] are simple 
matrix combinations of the matrix quotients 
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obtained from the continued fraction expansion or 
Routh's algorithm, 

2. The submatrices appearing below the main diagonal 
have the same value as the submatrices at the 
diagonal (with exception of [A 2 ]). 

3 . The submatrices which appear above the main diagonal 
can be expressed in terms of matrix quotients in a 
very regular way. 

4. State space formulations for the first and second 
Cauer form can be obtained directly from the Cauer 
mixed form through direct substitution. 

C. APPROXIMATION BY THREE CAUER FORMS 

The reduction of the order of a transfer function or 
decreasing the dimension of a state matrix is highly 
desirable or sometimes necessary in the analysis and design 
of control systems. 

In terms of continued fractions, the simplification 
problem is carried out by expanding a given transfer function 
into one of the Three Matrix Cauer Forms of continued frac- 
tion and ignoring some matrix quotients. If the given 
system is outlined in state variable form, the simplifica- 
tion method is realized by partitioning the matrix and 
discarding some parts. Several examples are included 
for demonstrating the power of the method. Also a thorough 
comparison among the three Cauer forms is presented and their 
advantages and disadvantages are discussed. 
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Feedback Gain and Feedforward Gain - Consider the system 
shown in Figure 5. The closed loop of the overall transfer 
function is known as follows: 

T(s) = [H 1 + [F ? + G]" 1 ]" 1 (42) 




Figure 5. Block diagram for a typical feedback system 
with a minor feedforward loop with two in- 
puts and two outputs. 

where G = 0. Equation (42) can be considered as the simplest 
continued matrix fraction expansion. The physical meaning 
implied in the formula is significant. It is easily seen 
that when is high the overall gain can be approximated by 
H" 1 , in other words, H-, dominates the behavior of the 
system. 
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This fact establishes the feedback loop as the most im- 
portant link for influencing the behavior of the system, 
leaving the feedforward loop as the second most important 
link. Furthermore if the subsystem G, instead of being zero 
is still a high order transfer function, it is possible to 
continue the expansion one after another. This corresponds 
to a combination of many feedback and feedforward blocks as 
shown in Figure 6. It should be noted that the most dominant 
term is and the second influence term is F^ . When the 
matrix quotients in the continued fraction are lower and 
lower in positions, they are less and less important as far 
as the influence to the performance of the system is con- 
cerned. This observation is the general basis for the sim- 
plification technique developed for multivariable systems. 

Considering a simple case, a second order transfer 
function such as: 



T(s) = 


[A 22 s + 


A 21"^ A 13 S + A 12 S + A ll^ 


(43) 


can be expanded into three different matrix Cauer forms 


of 


continued 


fraction 


as follows: 




F ., ( s ) 
cl 


= [H^S + 


• [H' 2 + [H’ 3 S + [H^]" 1 ]" 1 ]" 1 


(44) 


F ,(s) 

c2 


■ [H i + 


[H 2 | + CH 3 ♦ CH^]- 1 ]- 1 ]- 1 ]- 1 


(45) 


F As) 
c3 


* [K i + 


K'S + [K_- + Kl]" 1 ]" 1 
i IS 2 


(46) 
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where F , (s), F ~(s), and F 0 (s) are the first, second and 
cl c2 ’ c3 ’ 

mixed Cauer continued matrix fraction expansions of equation 

(43) . 

From previous observation, it is known that the most 
dominant term in equations (44), (45), and (46) is H^s , 
and + K'^s respectively. It is desirable to find a 
meaningful interpretation for the dominant term of each one 
of the expansions just performed. This task is accomplished 
by applying the initial and final value theorems to equations 

(44) , (45), and (46) as follows. 

1 . Cauer First Form 

Performing the inverse procedure on F c ^(s) it can be 
written as 



F d ( s ) =jj{fy= C H 2 H g H S + H 2 + H ^ ] [H^H^H^H^S 2 + (H^H^+H^H^+H^H^ ) S + I ] -1 



(47) 



Applying the final value theorem to equation (47) and 
allowing 



U(s) = 




1 

m 



where m is the number of inputs is 
found that 
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( 48 ) 



y(t-°o) = [H +H ] 
s+0 





similarly, applying the initial value theorem and allowing 



U( s ) 



1 

1 



1 

2 



1 

m 



y (t->0 ) = [H' ] _1 

S-J-oo 




(49) 



The meaning implied in equations (48) and (49) is very 
significant. The initial conditions dominate the behavior 
of the system. In other words, the Cauer first form 
influences very heavily the transient part of the response. 

2 . Cauer Second Form 

Performing the inverse procedure of ^^(s), this can 
be written as: 
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F c2 (s )=g[fy= CH 2 +H 4 S+H 2 H 3 H 4 ;ii:is2+(H l H 2 + H l H 4 + H 3 H 4 )S+H l H 2 H 3 H 4 :i 



-1 



(50) 

Applying the final value theorem to equation (50) and allowing 



U(s) 



m 



— it is found that y(t+«>) 
s 

s-*-0 



[H^] 



-1 



m 



(51) 



Similarly, applying the initial value theorem and allowing 







V 


1„ 




lo 


2 




2 


• 


y(t->0)=[H + H ] 


• 




s->°° 


# 


1 




1 


m_ 







(52) 



The results obtained in equations (51) and (52) imply that 
the final or steady state value dominates the behavior of 
the system. In other words, the Cauer second form influences 
very heavily the steady state part of the system response. 

3 . Cauer Mixed Form 

Performing the inverse procedure on F g(s) this can 
be written as: 

F c3 (s)= DTsT =[K 2 S + K 2 ][K, l K 2 s2+(I + K l K 2 + K i K 2 )S+K l K 2-l” 1 (53) 



34 



Applying the final and initial value theorems to equation 
(53) in the same fashion as for the previous forms, the 
following results are obtained: 



y ( t->°°) 
s-*-0 




and 



y ( t+0 ) 




1 



m 



(54) 



(55) 



These results show the steady state value and the initial 
conditions to dominate on equal levels of significance the 
behavior of the system. Thus the Cauer mixed form influences 
the system response in the transient part as well as in the 
steady state part. This fact makes the Cauer mixed form a 
better and more accurate device to be used in the simplifi- 
cation and reduction techniques of transfer functions of 
multivariable systems. 
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D. TRUNCATION OF THE CAUER MIXED FORM THE BEST APPROXIMATION 
IN A DEFINITE MATHEMATICAL SENSE FOR MULTIVARIABLE 
SYSTEM REDUCTION 

In a recent paper, M.J. Goldman and C.T. Leondes [3] 
developed the mathematical basis for the simplification 
technique involving the truncation of the mixed Cauer form 
for the single input single output case. In this section, 
an extension of their work to the multivariable case is 
presented and its superiority over any other form of con- 
tinued fractions is established. 

The transfer function matrix T(s) for multivariable 
system can be expressed as follows: 



T ( s ) = [A- S n 1 + A- .S n ' 2 +. . .+A_S 2 + A 00 S + A-, 3 x 
2 , n 2,n-l 23 22 21 



[A l,n + l Sn + VA"' 1 + -'- A 13 s2 + A 12 S + A li rl 



(56) 



where Ai , j are constant, m by m matrices and Aij = aj[I], 
j = l,2,...n+l, where each aj is a coefficient of the 
common-denominator polynomial or 

n+1 

A(s) = E ajS^ - and [I] is the identity matrix. 

3 = 1 

The nth convergent of a mixed matrix Cauer form can be 
represented by the following two configurations: 



36 



A B -1 = [K-, +KL S + S[K 0 +K' S+S[ . . .S[K +K' S] 1 1 1 ] X ] 1 
nnll l l nn 

A' B rl :-[K 1 -+K'+-[L-+K'+-[...-[K - + K* ] -1 ] -1 ] -1 ] -1 
nnslsls2s2s sns n 



(57) 

(58) 



where 



A = K A , + SK* A . + SA 0 , .. 
n n n-1 n n-1 n-2 o 



A =0 A ± = I 



B = K B + SK' B , + SB 0 , B =0 B, = K + KL S , 
n n n-1 n n-1 n-2’ o 111 



(59) 



A’ 



n 



-K A' , + K' A T n + -A _ , A' 
s n n-1 n n-1 s n-2 o 



= 0 A' 



ii 

s 



B’ = -K B’ , + K* B’ -i + -3' . , B' 
n s n n-1 n n-1 s n-2 o 



= I B' = 



— K + K' 

si 1 



(60) 



Recurrence relations (59) and (60) have been derived from 
standard results in the theory of continued fractions 
(Rice 19 64) and K^'s and K’^'s from the generalized matrix 
Routh Algorithm (Shieh and Gaudiano 1974). 

Since 



A B 1 = A' B' - 1 
n n n n 

then 



A = S n A' 
n n 



B = S n B’ n = 0,1,2, 
n n j j > 



(61) 



(62) 
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Proof : 



Substituting equation (62) into (59) 

S n A’ = S n-1 K A' , + SIC S n-1 A’ , + SS n_2 A 0 
n n n-1 n n-1 n-2 



and 



( 63 ) 



S n B’ = S n-1 K B’ . + SIC S n 1 B’ 0 

n n n-1 n n-2 

Dividing equation (63) by S n , recurrence relation (60) is 
obtained. (It is important to notice the generality of 
equations (59) and (60). 3y making all the K's = 0, the 
recurrence relationships for the nth convergents for the 
Cauer first form is obtained, similarly by making all the 
K'is = 0, the recurrence relationships for the nth covergents 
for the Cauer second form are obtained.) 



From (59) 



*-l 



and 



K = [A - SIC A , - SA 0 ]A . 

n n n n-1 n-2 n-1 



K = [B - SK' B , - SB 0 ]B 1 , 
n n n n-1 n-2 n-1 



Solving and simplifying terms 



(64) 



A B , - B A = -S ( A , B 0 - B , A „) 
n n-1 n n-1 n-1 n-2 n-1 n-2 



(65) 



so by induction 



A B , - B A = ( -S ) n-1 *I 
n n-1 n n-1 



where I is the Identity Matrix. 



( 66 ) 
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Hence, the difference between two consecutive convergents is 
given by 



A B 1 - A , B 1 ^ = ( -S ) n-1 [B v .B v . ,] -1 
n n n-1 n-1 K K-l 



(67) 



Similarly, from (60), and following the same procedure 



A’ B' , - B r A' , = ±(A' , B f 0 - B’ ..A' „ ) . 

n n-1 n n-1 s n-1 n-2 n-1 n-2 



( 68 ) 



So by induction 



B 

n n- 



1 



B’ A' 
n n-1 




n-1 



I. 



(69) 



and the difference between two consecutive convergents is 
given by 



A' B 71 - A’ Z 1 , = - * (- -) n 1 [ B’ B' 

n n n-1 n-1 s s n n-1 



-1 



(70) 



By looking at the recurrence relations (59) and (60), is 
observed that when S = 0 or S =o:> , respectively, all B^ and all 
B' are non-zero provided that all K and K' are non-zero. 
This is equivalent to apply the final value theorem to 
equation (67) for a unit step inputs; and applying the 
initial value theorem to equation (70) for an impulse 
response . 

Hence , 



d^ 

ds 



-[A B 1 
D n n 



A ,B -1 ,] C n 
n-1 n-1 S=0 



lim - I * S(-S) n 1 [3 B , ] 1 
s n n— -L 

s-»-0 

0 j= 0,1,..., n-2 

(71) 
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Similarly , 




lim S " - I (--) [B' B' , ] _1 
s s n n-1 



S-J-oo 



0 



j = 0,1 



n-2 



( 72 ) 



Equations (71) and (72) show the results that were expected, 
namely: The (n-l)th approximant and the nth approximant 

goes to zero in the sense of minimizing the follows semi- 
norms of the difference 



Note: A semi-norm is a norm which does not satisfy the norm 

axiom "||F(s) || = 0 implies F(s) = 0". 

Since the derivatives in (73) correspond to the coeffi- 
cients in the Taylor series expansion of the functions 

F ,(s) and F' ,(s) about the points S = 0 and S=°°. Then, it 
n-1 n-1 ^ 

can be deduced that the output difference to a nth approximant 




(73) 



where 



and 




(74) 
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has the same Taylor series as the original function for 

terms up to S 1 " 1- "*" and (i) n “‘'', effectively neglecting the part 

s 

of the transfer function which differentiates the input >_ n 
times in the steady state as well as in the transient. This 
gives the sense in which the approximation by the Cauer mixed 
form works. 

E. SIMPLIFYING A MATRIX TRANSFER FUNCTION 



If the 


following 


nth 


order system is 


given, 


T ( s ) = 




+ . 


. .a 24 s + a 23 s 


+ A 22 S + A 21 ] x 




CA l,n + l S " 


+ . 


..Ai 4 s 3 ♦ A 13 s 2 


+ A 12 d + A 1X ] 



(75) 

a simplified model of the system is desired. By performing 
the Generalized Routh's Matrix Array, the matrix quotients 
of the Cauer mixed form in equation (3) can be obtained. 

If an m order is desired only the first m pairs of K’s 
and K'is should be kept in equation (3) and the remaining 
should be omitted. After the inverse procedure has been 
performed on the truncated continuted fraction, the simpli- 
fied model is obtained. 

For example the general transfer function obtained by the 
matrix continued fraction expansion for the Cauer mixed form 
is given in equation (76). 

T(s) = [K 1 + K' 1 S+[K 2 |+K r 2 +[K 3 + K , 3 S+[K 4 |+K , 4 +[ . . . I" 1 ]" 1 ] -1 ]" 1 ]" 1 
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If a second order simplified model is required, only the 
first two pairs of matrix quotients will be kept, that is 



K 



, k ' i5 



1 < 2 » and K' 2 , and the rest should be discarded. 



After the inverse procedure has been performed on the 
truncated matrix continued fraction expansion, the transfer 
function given in equation (77) is obtained. Equation (77) 
is the simplified model in Mixed Matrix Cauer Form. 



T ( s ) = [K' 2 S + K 2 ] x [K , 1 K , 2 S 2 +(I+K 1 K' 2 +K , 1 K 2 )S + K 1 K 2 ]" 1 



(77) 



The truncation procedure outlined above, applies in similar 
manner for the Cauer' s first and Cauer ' s second forms. 

This methodology is particularly advantageous when state 
space terminology is used. For the Cauer mixed form, the 
state space formulation is written as follows: 



KI K' 


Ki K' 




Ki 1C, 






. . . Ki K' 




1 2 


1 4 




1 6 






1 n 




K 1 K 2 


(K i tK 3 


)K, 4 


(k' 1+ k' 3 


)K' 

6 




. . . (K^ + K'g 


)K' 

n 


KI K' 


(k;+k; 


)KT. 


(K'+K' 


+ K'r 


) K'_ 


. . . (Ki +K'_ 


+K')K' 


1 2 


1 3 


4 


1 3 


5 


6 


1 3 


5 n 


_ K i K 2 


(K i tK 3 


,K 4 


(K'pK'3 


+K 's 


)K* 

6 


. . . (k'^+k' 3 


+KI+ . . .+K' 
5 n 


I + ^K 


2 +K l K 2 


K^K 


4 +K l K 4 






. . K' x 


K +K, K' 
n In 



k , 1 k 2 +k 1 k' 2 

K 1 K 2 +K 1 K 2 



K i K 2 +K l K 2 



I+(K^+K' 3 )K 4 +(K 1 +K 3 )K' 4 
(K , 1 +K , 3 )K 4 + (K 1 + K 3 )K[ | 



(K! l +K , 3 )K 4 +(K 1 + K 3 )K' 4 



n-1 n 



(K'+K' )K + ( K, +K, ) K' 
1 3 n 1 3 n 

(K» +Kh+KV)K + 

1 J 5 n 

(Kt+K-+K c )K 
. 1 3 5 n 



«i +K ’ 3 + 



+ K' . )K + 
n-1 n 



(K, +K-+ . . . +K , ) K' 
1 3 n-1 n 



n/2_ 

E, 



J n/2 
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K 1 K 2 K 1 K 4 ••• K l K n 




E 1 




I 


K 1 K 2 (K 1 +K 3 )K 6 ••■«l +K 3 ,K n 




E 2 




I 


K 1 K 2 (K 1 +K 3 )K 4 <K 1 + W K 6 ' ' ‘ <K l + W K n 

• • • * 




E 3 


+ 


I 


• • • * 

_ K 1 K 2 <K 1 +K 3 ,K 4 ( W K 5 ,K 6 ••• <K l +K 3 + -" +K n-l )K n 




fn/2_ 




I 



(78) 



and 











E i 




E 2 






E 2 


Y = [K 2 , K 4 , K 6> ... K n 3 


E 3 


+ [K’ 2 , K’ 4 , Kg , ... 


K' ] 
n 


E 3 




_ E n/2_ 






_ E n/ 2_ 



(79) 



where 



6 1 




^m+1 




e (n/2-l)m+l 

f • 


S 2 




e m+2 


# , 


e (n/2-l)m+2 


• 

• • 

e 

_ mj 


E 2 * 


_ e 2m_. 


E n/2 


• 

♦ * 

_f (n/ 2 )m 



I = identity matrix, and E^, . .. £ n /2 > E 2 > ^n/2 

have the same dimensions as E^, E^, . .. E n ^ 2 respectively 
and m corresponds to number of inputs or outputs. 
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It is interesting to note that the simplification of the 
state equations can be carried out by partitioning the matrix, 
or by only keeping a part of the original matrix as a simpli- 
fied model. In other words, if a simplified model of a two 
dimensional matrix is required which is equivalent to a 
2 x (2m) order transfer function, the upper left hand corner 
of the original matrix is taken as the simplified model. 



Therefore : 



K'K' K'K’ 




e„ 




I+K! 


, K 0 +K, KL K’K.+K-.Kl 




• 

E, 


12 14 




l 


- _ 




L 2 1 2 14 14 




1 


K 1 K 2 (K i +K 3 )K 4 




r 

_^2_ 




_ K i K : 


? + K 1 K*2 I+(K , 1 +K , 3 )K 4 +(K 1 +K 3 )K' 4 






K,K 0 K..K, 




E, 




i 






12 14 




1 


+ 




u 




K 1 K 2 (K 1 +K 3 )K 4 








i 


(30) 



Y = CK 2 , K 4 ] 


e i 

E n 


+ nc 2 , K’ 4 ] 


E 1 

e: 




2 




2 



where the vector E^, , E^ , E^, E^, E 2 , keeping the same 

dimensions as were given for the equations (78) and (79). 

The method used for the simplifications of a transfer func- 
tion and state space formulation based on the Cauer mixed 
form holds in the same fashion for the Cauer first and second 
forms . 

For effects of simulation Figure 7 shows a diagram for 
two inputs two outputs fourth order system for the Cauer 
Matrix mixed form, in similar form for the Cauer first and 
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Figure 7. Simulation diagram for fourth order system with 
for Cauer matrix mixed form. 



m * 2 
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second form a simulation diagram can be obtained as given 
in Figures 8 and 9. 

It is important to notice the number of reductions, q, 
which one can realize are constrained by the number of inputs, 
m, and the order of the system, n, where q = n-km for 
k = 1,2,... and q > 0 . 

E. NUMERICAL EXAMPLES AND GRAPHICAL COMPARISON OF THE 
THREE CAUER FORMS 

1 . Example 2 . 1 

Consider two inputs, two outputs fourth order system 
where a reduced second order system is required. 

The state space equations for the system are: 



The 





0 


0 


3 .25 




L.125 




0 


0 




0 


0 


2.25 


-( 


3.125 




0 


0 


• 

X = 












x + 








-1 


0 . 5 


-1.5 


-0 . 5 




0 . 5 


0 




-1.5 


-0 . 5 


-1.0 




+ . 0 




1 


0 . 5 




1 


-1 


1 


1 










Y = 


•0 


1 


-1 


1 


X 









transfer function for this system 



is given by equation 



(82) 



(83) 
(84) . 



. — - 




— 


s + 1 s — 1 




2s 2 + 3s + 2 s — 1 


-2 s + 1 




-4s 2 - 4s - 1 2s 2 + 6s + 3 



( 84 ) 
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Figure 8. Simulation diagram for fourth order system with m = 2 
for Cauer matrix first form. 



4 7 




Figure 9. Simulation diagram for fourth order system with m = 2 
for Cauer matrix second form. 
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From the transfer function, the matrix Routh's Array and 
the following matrix quotients for the three Cauer forms are 
obtained . 

Therefore, reduced second order models are: 

Cauer first form 




H’ 



5 

1 



-1 

5 



^L_ 

13 



Cauer second form 





2 1 




0.5 0 


r 




H 0 = 






1 4_ 


2 


_0 . 125 0 . 2 5_ 



Cauer mixed form 




where their respective transfer functions are given in equa- 
tions (84), (85), and (86). 

Cauer first form: 



F cl (s) ■ 



3c + 2.5 
s + 0 . 5 



Cauer second form: 



s - 0 . 5 
s + 2.5 



2s 2 + 10s + 6.5 



(84) 



F c 2 (s) = 



4s + 2 
s + 1 



-1 

2s + 2 



9s + 13s + 5 



(85) 
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Cauer mixed form 



F c3 (s) 



3s + 2 
s + 1 



s - 1 
s + 2 



1 

2s 2 + 3s + 5 



( 86 ) 



Figure 10 presents the step responses of reduced models just 
developed, where (10a) corresponds to the output (1) and (10b) 
corresponds to the output (2). 

2 . Example 2 . 2 

The following example shows the reduced systems in 
Cauer first and second forms are only stable for a second 
order approximation, their fourth order approximation pro- 
vides an unstable response. Whereas, the fourth order 
approximation for the Cauer mixed form is completely stable. 
Furthermore, the Cauer first form gives a poor approximation 
in the steady state portion as is expected. Consider the 
following sixth order system with two inputs two outputs as 
given in equations (87) and (88). 



0 


0 


1.0 


0 


0 


0 




0 


0 


0 


0 


0 


1.0 


0 


0 




0 


0 


0 


0 


0 


0 


1.0 


0 


x + 


0 


0 


0 


0 ' 


0 


0 


0 


1.0 




0 


0 


-50 


9.96 


-108.5 


19 . 6 


-16.5 


4 




2 


-1 


49 


-19.96 


107 


-39.6 


13 


-10 




-2 


2 



(87) 



10 



0.4 1 0.2 0 0.2 



Y 



x . 



1.0 4 21 



1 0 



( 88 ) 
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Figure 10. Outputs for the original fourth order system and reduced second order system in the three 
Cauer forms. 



The transfer function of this system is given by equation (89). 



T ( s ) = 



s + 10 



0.2s + 0.2s + 0.1 



S +2s+l s+4 



s 3 + 10s 2 + 55s + 25 



0.5s 2 + s 2 + 0.2s +0.2 



s 3 + 3.5s 2 + 1.5s + 0.5 s 3 + 6s 2 +20 s + 10 



-1 



or if the inverse is to performed 



(89) 



T ( s ) = 



_ -0.25s 5 + 0.1s 4 + 14.9s 3 + 79.25s 2 + 209.75s + 99.95 
_s 5 + 7s 4 + 25.5s 3 + 40.5s 2 + 33.5s + 8 
0.2s 5 + 1.7s 4 + 7.1s 3 + 6.8s 2 + 8.48s + 2.3 
-0.5s 5 - s 4 +11. 3s 3 + 93 . 58s 2 + 244.76s + 99.98 



0.5s 6 + 13.25s 5 + 130.55s 4 + 562.53s 3 + 1349 . 13s 2 +1049 . 87s+249 . 99 



(90) 



Cauer Matrix First Form 

Performing the respective Routh's Array: 



ri 0.5 ~i rio 

|_1 1 J |_3 . 5 


1 1 

i — 1 CD 


ro 0.2 in 

[_i o J |_2 


0.2 J 


["5.5 -0.5 1f29 
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{_- 3.5 4 J |_-49.5 


15 . 5 


(“2.686 -0.5 


9.1 


-0.444] 


[_-2.5 1.417 


_-3 . 95 
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“2 5 


0.02 "I 


[-53.664 3 .163 


_0. 5 


10 J 



25 

0.5 



0 . 2 
20 

0.1 

4 

0 . 02 " 
10 
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7.234 -0.055 

3.046 3.696_ 

"2 5 0.02 “ 

0.5 10 



and from the Routh's array, the respective H’i's are: 





2.5 1 




0.035 0.054 






H' = 








2 






_5 1_ 




_0 . 19 7 0 . 0 2 5_ 



2 .566 


0.556 


H! = 


0.076 


-0 .039 






4 






_1 . 9 8 6 


3 . 5 26_ 




_-0. 28 


0.006 



1.64 


-1.327 


H' = 


0 . 29 


-0.006 






6 






_-7.7 3 


0.74 




_0 . 114 


0 . 3 6 9 _ 



The state space representation in Cauer first 
is given by equations (91) and (92). 
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0 




0 1 








0 0 


0 


0 


1 


0 




1 0 
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0 1 
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y 



0.035 0.054 

0.197 0.024 



0.076 -0.039 | 0.289 -0.006 

I z 

-0.28 0.006 | 0.114 0 . 3 6 9_ 

(92) 



The reduced fourth and second order model are obtained by 
partition of the respective matrices in equations (91) and 
(92). For this case, the fourth order model is unstable and 
the second order model gives an untolerable error in the 
steady state portion of the response, so this method is not 
applicable for this example . 

Figure 11 shows the step response for the original 
system and the reduced second order model, 
b. Cauer Matrix Second Form 

Performing the respective Routh's Array: 
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0.02 


f55 
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0 . 5 
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Figure 11. Outputs for the original sixth order system and reduced second order system by Cauer first form. 



and from the Routh's array the respective Hi's are: 





2.505 


-0.058 




0.191 


0.003 


H 1 5 


_-0 .2 


2 . 5 0 5 _ 


H 2 * 


_0 . 03 3 


0 . 2 2 9_ 




”-71.026 


-9.128” 




”0.013 


-0.038 - 


H 3 * 


_- 8 2 . 5 2 7 


-55 .455_ 


H 4 * 


_0 . 083 


0 . 06 5_ 




" 92.956 


25.02 ” 




”0.604 


0 .429" 


H 5 * 


145.353 


67 .357 


H 6 * 


2 .05 


-1.293 



The state space representation in Cauer second form is 
given by equations (93) and (94). 



-0 .476 
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-3 .263 


0.322 


8 . 3 


1 

LD 

CM 

OO 

1 




0 


1 



Y 



0.191 0.008 

0.033 0.229 



0.013 -0.038 | -0.6 0 .‘429 

I 
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(93) 



(94) 



The reduced fourth and second order model are 
obtained by partition of the respective matrices in equations 
(93) and (94). For this case like the Cauer first form, the 
fourth order model is unstable and the second order model 
offers a good approximation to the step response of the 
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original system. Figure 12 shows the outputs of the original 
and reduced second order model for a step input, 
c . Cauer Matrix Mixed Form 

Performing the respective Routh's Array: 



25 0. 

_0 . 5 10 


02] p55 

J Li-5 


o. 2 i rio 
20 J L 3 * 


"10 0 . 

_1 4 
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-4.493 “| 


“5.55 8 
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and from the Routh's array the 




.04 



respective Ki's and K'i's are: 
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K' 2 = 
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. 061_ 
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-1, 


. 365“ 


K 3 
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K, 3 * 
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2 , 
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The transfer function for the reduced fourth order system 
is given in equation (95). For this model the step response 
is completely stable and the approximation has an error of 
only 0.5% over the entire response of the system. Figure 13 
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Figure 12. Outputs for the original sixth and reduced second order system by Cauer second form. 
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Figure 13. Outputs for the original sixth and reduced fourth order system by Cauer mixed form. 



shows the outputs of the original and reduced fourth order 
model for a step input. 



F(s) 



-0.012s 3 - 0.063s 2 + 1.44s + 0.972 
_ 0.061s 3 + 0.773s 2 + 0.278s + 0.078 
0.012s 3 + 0.173s 2 + 0.068s + 0.022“ 
-0.03s 3 - 0.234s 2 + 1.784s + 0.972 



1 

0.03s 4 + 0.9795s 3 + 8.215s 2 + 8.718s + 2.431 



(95) 
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III. DERIVATION OF LOWER ORDER SYSTEM FOR THE LINEAR 
REGULATOR PROBLEM VIA CAUER FORM 

A. LINEAR REGULATOR SYSTEM 

R. E. Kalman [12] has shown that when a system is 
described by equations ( 96 ) with a performance function 
given by equation ( 97 ) 

x = Ax + Bu ( 96 ) 

J = J 5E[x^Qx + u^Ru]dt (97) 

and if the plant is completely controllable then an optimal 
u" exists and is given by equation ( 98 ) 

u* = -Kx ( 98 ) 

where the feedback matrix X for the controller depicted in 
Figure 14 is a constant matrix as t-*». 

In recent papers, Goldman [10] and Aoki [18] have 
derived a way to obtain near optimal solutions using the 
"aggregation matrix" , to transfer an arbitrary state x to the 
origin of the state space while minimizing the criterion 
function given in equation ( 97 ) . 

The present work will show that the reduced order optimal 
regulator can be obtained where the original system is 
translated in Cauer form and also a near optimal solution 
can be found through a mere partition of its optimal solution. 
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PLANT 




Figure 14 . Linear Regulator System 



B. THE STATE VARIABLE REPRESENTATION BY CAUER 
SECOND FORM AND ITS OPTIMAL SOLUTION 

It is known that the standard state variable representa- 
tion of a linear time-invariant system is given by equations 
( 99 ) and (100 ) . 

x = Ax + Bu ( 99 ) 

y = Cx (100) 



T 

where x is a column vector nxl given by [x n , x„. . .x ] , A 

J 1 2 n 

is an nxn constant coefficients matrix, x is a column vector 

T 

nxl given by [x 1 , x 2’’' ,X n^ , B is an nxl constant coefficients 
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matrix, u is the input, y is the output and C is lxn constant 
coefficients matrix. 

Goldman (10), Chin and Shieh (11) have proven that the 
system given by equations (99 ) and (100) can be represented 
in Cauer second form as shown in equation (101) and (102). 



i - Fz + Gu 


(101 ) 


y = Mz 


(102 ) 


where , 




F = PAP -1 


(103 ) 


G = PB 


(104 ) 


M = CP" 1 


(105 ) 



where the matrix P is an nxn upper triangular matrix and the 
elements in the triangle are copied directly from the elements 
of the Routh’s array, where the nth row of the P matrix 
is the (2n+l)th row of the Routh’s array and the states 
variables x and z are related by equation (106). 

z = Px (106 ) 

If the system in ( 99 ) and (100 ) is put in transfer function 
notation the result is given by (107 ) • 

m / s_a21+a22S+a23S 2 +a24S 3 + . . .+a2nS n_1 , nn _s 

T(s) * s — (107) 

all+al2S+al3S +al4S +. . .+a2n+lS 
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The Routh's array can be formed as 



all 


al2 a!3 


. . . aln al ,n+l 




a21 


a2 2 a2 3 


. . . a2n 


( 108) 



and the elements of third, fourth, and subsequent rows can 
be evaluated from the following algorithm. 



ajk = aj-2, k+l-f j -2aj -1 , k+1 j=3 ,4 , . . . ,n+l 

k = l , 2 , . . . 



fp 



ap,l 
ap+1 , 1 



P = l,2 , • . -n 



(10 9 ) 
(110 ) 



The elements of the matrix F can be obtained also from 
equation (111) where ap,l and ap+1,1 correspond to the 
elements of the first column of the Routh's array and then 
F is formed as shown in equation (111). 





f 2fl 


f 4f 1 


f 6f 1 


. . f nf 1 




f 2f 1 


f 4 ( f 1 + f 3 ) 


f 6 ( f 1+f 3 ) 


. .fn(fl+f3 ) 


F = - 


f 2 f 1 


f 4 ( f 1+f 3 ) 


f 6 (f 1+f 3+f 5 ) . 


. . fn( fl+f 3+f 5 ) 




_f 2f 1 


f 4 ( f 1 + f 3 ) 


f 6 ( fl+f 3+f 5 ) . 


. .fn(fl+f3+. . .fn-1) 



(Ill) 



The elements of the matrix G correspond to a column vector 
of nxl of l's and the elements of the matrix M can be formed 
also by fp's where p=2 ,4 , 6 , . . . 2n as shown in equation (112). 

M = [ f 2 , f 4 , f6 , . . . f 2n] 

(112 ) 
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The performance function given in equation ( 97 ) can be 
translated to Cauer form by using the relationship given 



in equation (106), i.e., 

J = VCz* T Q 2 “ + u* T RU*]dt 
c c 

where 

Q = P" 1T Q P _1 
c x 

and the optimal u* is given by equation (115 ) . 

u* = -K z* 
c 

where 

K = K P " 1 
c 

So the equivalent system in Cauer form given in Figure 14 

can be represented by Figure 15 . 




Figure 15 . Linear Regulator in Cauer Form. 



(113 ) 
(114 ) 

(115 ) 
(116) 
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Note that K is the feedback gain matrix of the Cauer 
c 

system and u* is the optimal control law of both systems 
(Figure 14 and 15). 

From the theory of optimal control, z* and u* are found 
by solving the set of necessary conditions: 

z*= (z*, u* , p*, t) 

3p* 



P :: = ~ 



dH 

3z ; * 



( z" 



U‘ 



t) 



0 = 



dH 



du 



■x- ( Z * , U*, p* , t) 



for all 
te [f 0 , f f ] 



n _ d 2 H . 

0 ‘ (z "’ u “’ p “’ 



t) 



(117 ) 



Where the function H is called the Hamiltonion and is 
defined by: 



tf(z,u,p,t)=%[z T Q c Z+U T RU]+p T [Fz+Gu] 

(118 ) 

and p is the Lagrange multiplier or Costate state. 

Substituting equation (118 ) into equation (117 ) , the 
canonical system in equation (119) is formed. 





G 



T 




(119 ) 
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Fortunately, for the Optimal Regulator problem, it is not 
necessary to solve these equations. Kalman [12], Tyler, 

J.S. and Tuteur, F.B. [19] have shown that when the optimal 
control u* is generated by equation (115), the solution 
for is obtained from equation (120) where the matrix, T, 
is the solution to equation (121) in steady state. 

K c = R~ 1 G T T (120) 

0 = f t t+tf-tgr _1 g t t+q c (121) 

Equation (121) is the steady state form of the Ricatti 
equation . 

For the single input case the weighting matrix, R, is 
a scalar. When this is the case, a lower order optimal 
linear regulator can be found just by mere partition of 
F, G, M, and K c matrices as shown in equation (122). 

Hence the new system is depicted in Figure 16 . 




Figure 16 . Reduced Linear Regulator in Cauer Form 
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H cnn 



(12 2 ) 

The reduction scheme outline above has left u * and u 

c r 

unchanged. Therefore, u p " is the suboptimal solution for 
all three systems, that is the original system, the trans- 
formed system and the reduced system. 

C. THE REDUCED ORDER SOLUTION 

Based on equation (122), the reduced system is given by 

zr = Frz + Gr u (123) 

y = Mrz (124) 

where Fr, Gr and Mr are the partitions of the matrices F, G 
and M respectively. 
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Goldman [10] has shown that the original system is 
related to the reduced order model by the following equations: 



zr = eAe + zr + eBu (125 ) 

y = Ce + zr (126 ) 

x = e + zr (127) 

zr = ex (128 ) 



where the matrices e and e are partition of the P and P 
matrices in rectangular form respectively, such that e is of 
order (rxn) and e + is (nxr) , where, n, is the order of the 
original system and, r, is the order of the reduced system. 

From equation (127 ) and ( 98 ) it is easy to show that 

u* = Ke + zr (12 9 ) 

then 

Kr = Ke + (130 ) 

or 

K = Kre (131 ) 

Note that from equations (98 ) and (106 ) 

K = KP -1 (132) 

c 

or 

K = K P (133) 
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Then from equations (130) and (133) 

Kr = K Ir (134) 

c 

where 

Ir = Pe + (135) 



The matrix Ir has dimension (nxr) and is a special type of 
identity matrix where the last (n-r) rows are zeros and the 
first rxr components correspond to the identity matrix, in 
the same way, it is possible to show that the different 
relationships that exist for the weighting matrices, Q, Qc 
and Qr are given in equations (136), (137), (138) and (139). 



-IT -1 
Qc = P Q P 



+T~ + 

Qr = e Qe 



+T T + 
Qr = e P QcPs 



Qr = Ir Qc Ir 



(136) 

(137) 

(138) 

(139) 



1 . Example 3 . 1 

The transfer function of a system is given in equation 
(]_40). It is desired to control the system in such a way that 
a performance function given in equation (141) is minimized. 

Due to the systems complexity, a near optimal solution is 
desired . 

H(s) = — i (140) 

S J + 6 .IS +5. 6S + 0 . 5 
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J = A[x x Qx+U T RU]dt 



(141 ) 



where 



5 0 0 



Q = 



0 4 0 



0 0 



1 



R = 1 



From the transfer function 
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CO 
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C= [1 0 0] 



The Routh ' s array 
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From the Routh' s array 





5.6 


6.1 


1.0 


p = 


0 


5 . 18 2 


1.0 




0 


0 


ol 

I — 1 



(142 ) 
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then 
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0 
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The Qc matrix obtained from equation (136) is 



”0.159 


-0.188 | 
i 


0 .029 


-0.188 


1 

0.372 j 
1 


-0.183 


0 .029 


-0.183 


1.154 



(143 ) 



(144) 



(145) 



(146) 



The Fc, Gc and Me matrices from equations (103), (104) and 
(105) are 





-0 .089 


0.105 


i 

i 

I 


0 .015 






Fc = 


-0.089 


-0.978 


i 

i 

i 

i 


0.139 








_-0 . 0 8 9 


-0.978 




-5 . 031_ 




(147) 


Gc = 


”1“ 
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0 

1 1 

II 

a 
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1786 


-0 . 211 


0 .032] 





1 
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The reduced system is given by equation (148) and the 
corresponding matrices that describe the system in equation 
(149) . 

TTll/ N -0 .0316s + 0 .193 

H4(s) = — Tj 

s +1.065s+0 .0965 (148) 



-0 . 089 



Fr= 



-0.089 



0 .105 
-0.978 



1 



Gr = 



1 



Mr= [0.1786-0.211] 



0.159 -0.188 



Qr = 



-0.188 0.372 



(149 ) 



The computer solution for the three systems for the feedback 
gain are shown in Table I. 



TABLE I 



SYSTEMS 


OPTIMAL FEEDBACK MATRIX K 


ORIGINAL 


1.79128425 


2 .087266 


0.41034 


CAUER FORM 


0 . 319649369 


0.0243818 


0.06619982 


REDUCED CAUER FORM 


0.31965241 


0 .02258695 
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Table II shows values found from the reduced system 
back to Cauer form and original system with the equation 
(131) and (132). 



TABLE II 



SYSTEMS 


NEAR OPTIMAL FEEDBACK MATRIX 


ORIGINAL 


1.7908 


2 .067 


0 . 342 


CAUER FORM 


0.31965 


0 .022586 


0.0 


REDUCED CAUER FORM 


0 . 31965 


0.022586 





Note that a near optimal feedback matrix is obtained directly 
from equation (131) but not to the Cauer system since the 
reduction technique discarded the last terms . 

Figure 17 shows the step response for the optimal and 
reduced linear regulator. Figure 13 shows the optimal and 
suboptimal control laws. Figure 19 shows the step response 
for the optimal and suboptimal linear regulator. 

2 . Example 3 . 2 

The methodology just developed is applied to a simple 
model of a nuclear reactor power generator [20, 21]. 

The heat generating process of a nuclear reactor is 
dependent upon the mechanism called fission (a fragmentation 
of matter) . The power generated by this process is directly 
related to the population of neutrons, n(t) and can be 
described by the following differential equation (developed 
from a diffusion balance equation) . 
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FIG. 17 - EXAMPLE 3.J. STEP RESPONSE 

OPTIMAL AND REDUCED LINEAR REGULATOR 




XSCflLE= 2.00 UNITS/INCH 

YSCflLE= 0.04 UNITS/ TNCH 



RUN NO. 1 
PLOT NO. 1 
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FIG. is - EXAMPLE 3,1 - NEAR OPT. SOLUTION 
OPTIMAL AND SUBOPTIMAL CONTRL LAWS 




XSCALE= 2.00 UNITS/INCH RUN NO. 1 

YSCALE= 0.20 UNITS/INCH PLOT NO. 1 
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FIG. 19 - EXAMPLE 3.1 STEP RESPONSE 

OPTIMAL AND SUBOPTIMAL LINEAR REGULATOR 




X5CALE = 2.00 UNITS/INCH RUN NO. 1 

YSCALE= 0.04 UNITS/INCH PLOT NO. 1 
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(150 ) 



n(t) = (IMtM 



) n(t) + Ac(t) 



c(t) = (j)n(t) - Ac(t) 



( 151 ) 



where 



6k(t) = 5kc(t)-an(t) 



(152 ) 



The variable 5k(t) is the input to the process and is given 
the name "reactivity" . It is clear by inspection that 
$k(t)<3 for stable system (in a linear sense). The 
variable c(t) is a measure of the concentration of fragments 
(precursors) that produce delayed neutrons according to a 
time delay (1/A) called the "half-life" of the precursor. 

The input, 5kc(t), is the control (reactivity) that is 
associated with the control rod position and, a, is a 
temperature feedback (reactivity) coefficient. 



at t=0 n(0)=10kw (as the operating output in steady 

state conditions) 

The control problem is stated as : 

Find the Optimal Control Policy u*(t) that will 
transfer the power level n(t) from the operating level 



The parameters for this system are: 



a = 10 - ^ kw 

_3 

1 = 10 sec. 



6 = 0.0065 



A = 0.1 sec 



-1 
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n(0) = lOkw to <3. new level n(t) = 50 kw where 
u( t ) = 5 kc ( t ) 



( 153 ) 



using the performance measure given in equation ( 97 ) where 



10 



-6 



Q 



0 



0 0 
0 0 
0 25 



(154) 



R = 1 



end compare this solution to the near optimal one obtained 
from the reduced order model. 

From the state of the problem the following diagram 

is drawn 




Figure 20. Block Diagram Optimal Control of a Nuclear 
Reactor . 
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The double subscript, ss, means steady state and the 
new definition of state variables are: 



xl(t) 


= n(t) - 


n 


(155) 


ss 




x2 ( t) 


= c ( t ) - 


c 

ss 


(156 ) 


x3 ( t ) 


= <Skc(t) 


- Skc 


(157) 



ss 



The steady state conditions are 
n(t) = 0 



(158) 



c(t) = 0 (159) 

<5 k( t ) = 0 (160) 

The initial conditions at n(t) = 50kw are: 

from the state variables definitions, equations (155), 

( 1 5 6 ) and (157) and equations (150), (151) and (152). 

c(t) = — n(t) - c(t) 



c( 0) 
c( 0 ) 



it n(0) 



640 



(161) 

(162) 

(163) 



x2 ( 0 ) 


=c(0) - c 

ss 




(164) 


x2(0) 


= -2560 




(165) 


5 k(t) 


= an(t) - XI 


°<t). + s 
HTtT 3 


(166) 


6 k( 0 ) 


II 

P 

2 

0 

v_/ 

1 


C(0) n 

nToT 3 


(167) 
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6 k( 0 ) 



(168 ) 



10 



-4 



x3 ( 0 ) 
x3 ( 0 ) 
xl ( 0 ) 
xl(0) 



dn 



ss 



XI 



c 

ss 



n 



ss 



+ 3 



5x10 



<S k( 0 ) -5k 

ss 

-4 

-4x10 
n( 0 ) -n 

ss 

-40 



then 





xkot 




-40 


x( 0 ) = 


x2 ( 0 ) 


= 


-2560 




_x3 ( 0 )_ 




— 4 

j-4xl0 _ 



and the optimal control law U*(t) will be 
u*(t) = -Kx ( t ) 



(169 ) 
(170 ) 
(171) 
(172 ) 
(173) 
(174 ) 



(175) 



same as given in equation ( 98 ). The system in variable 
form is given by equations ( 99 ) and (100). Since equations 
(150) through (157) are non-linear linearization is required 
(see appendix A). The linearized system is represented by 
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- i( an +6 ) 
1 ss 


X 


n 

ss 

i 




0 


x(t) = 


s /i 


-X 


0 


x + 


0 




0 


0 


0 _ 




_ 1 _ 



u 



(176 ) 



y(t) 

where 

A = 



= [1 0 0]x( t ) 



-6.9 0.1 5x10 

6.4 -0.1 0 

0 0 0 



( 177) 



( 178) 



B = 



0 

0 

1 



(179) 



C = [1 0 0] 



(130) 



Solving the Ricatti equations, the matrix K is given by 
equation ( 181 ) . 



K = [2.5x10 4 8xl0 -6 7.05] 



(131) 



The transfer function for the optimal system will be: 



H(s) 



5x10 4 (S+0.1) 

S 3 +14.05S 2 +61.9S+4.16 



(182) 
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The Routh's array for this system is: 



4.16 61.9 14.05 1.0 




From the Routh's array 



"20.3 


14.05 


l.cf 






0.0 


14.157 


1.0 






_ 0 . 0 


0.0 


1 . 0 _ 




(183 ) 


"0.049 


-0.049 


-0.0004 




0 . 0 


0.07 


r- 

o 

o 

1 




0 . 0 


0 . 0 


1.0 




(184 ) 



and the rectangular partition of P and P matrices give 



20 . 3 


14.05 


1.0 


_ 0,0 


14.157 


1 . 0_ 


“0.04 9 


-0 . 04 9" 




0.0 


0.07 




0.0 


o 

o 





(185) 



(186 ) 



83 



and from equations (125) and (125) the reduced second order 
suboptimal system will be 



« 


-0.205 


-2.735 




1 


zr(t) = 


-0 .205 


-4.169 


zr ( t ) + 


1 



y ( t ) = [246.305 3287 . 5]zr( t) 

The corresponding transfer function is given in equation 
(189 > • 

Hrop ( s ) 

S +4 . 3738S+0 . 294 

In Figure 21 is shown the step response for the original 
optimal and reduced suboptimal system. 



(187 ) 
(188 ) 

(139 ) 
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00 4.00 8.00 12.00 16.00 20.00 24.00 28.00 32.00 36.00 40.00 44.00 48.00 



FIG. 21 - EXAMPLE 3.2 _ STEP RESPONSE NUC. REP.C. 

ORIGINRL RND REDUCED OPTIMAL SYSTEMS 




T.QQ 



.00 



1S.OO 



24.00 



32.00 



40.00 



48.00 



56.00 



.00 



” 72.00 



XSCALE= 8.00 UNITS/ [ N C H 

TSCALE= 4. 00 UNITS/INCH 



RUN NO. 1 
PLOT NO. 1 



85 



IV. CONCLUSIONS AND RECOMMENDATIONS 



Comparative analysis of Cauer forms methodology for 
multivariable system reduction is established. The developed 
methodology is based on the Cauer matrix generalized form, 
which offers the closest approximation to the original 
system. The proposed methodology and because of the nature 
of the Cauer mixed form shows it to be superior to any other 
method proposed to date since it provides satisfactory 
results for both the transient and the steady state portion 
of the system response. The methodologies in state space 
as well as in the S domain are developed. For the basis 
of comparison, reduced order models using the three Cauers 
forms are obtained for two different examples . The results 
clearly show the superiority of the Cauer mixed form over 
the entire frequency range of system's response. The 
proposed methodology is algorithmic therefore, it is 
amenable to digital computation. 

A lower order optimal linear regulator can be obtained 
by mere translation of the original system to Cauer second 
form and their partition of the different matrices as shown 
in Section III-B. A suboptimal feedback matrix for the 
original system can be obtained by multiplication of the 
lower feedback optimal matrix by the rectangular partition 



86 



of the matrix P. The responses of both systems, the original 
and the reduced, as well as optimal and sub-optimal systems 
are extremely close. 

The results presented here are encouraging, there is 
a distinct need for future research, in particular in the 
reduction of systems with any given number of inputs or 
outputs . 



87 



APPENDIX A 



Linearization of Multivariable Systems 

Let equations (A-l) and (A-2) be the state variable 
modeling 



x = F(x,u) 



(A-l) 



y = g(x,u) 



(A-2) 



For the system described by equations (A-l) and (A-2) is 

operating at steady state conditions ( constant input u , 

s s 

producing constant state x and constant output y ) . The 

s s s s 

combination of these produces 



0 = f (x , u ) 
ss ss 



y = g(x , u ) 
^ss ss ss 



( A- 3 ) 
( A-4 ) 



If the system is perturbed by either drawing the states or 
the inputs, the system motion satisfies 



6 x = F ( x t5x, u +<5u) 



ss 



ss 



y t5 y = g ( x +5 x , u +6u) 
^ss 7 6 ss ’ ss 



(A-5 ) 
(A- 6 ) 



Both functions, F and g can be expanded in a Taylor series 

about the points (x , u ) resulting in the following 

s s s s 

representation of the system equations. 
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x + 8x=F(x , u ) +AS x-^5 u+a(5 x ,5 u) 
ss ss’ ss 

y +6y=g(x , u ) +G5 x+D5 u+8 (5 x ,5 u) 

ss ° ss ss 

where : 

A = nxn matrix 
B = nxm matrix 
C = vxn matrix 
D = vxm matrix 

n = order of the system, m = number on inputs, v=number of 
outputs. 

The functions a(5x,du) and S(5x,6u) represents all 

second order and higher terms in the Taylor series expansion. 

Substitution of equations (A-l) and (A-2) evaluated at 

(x , u ) and neglecting second order terms and higher 
s s s s 

terms, yields the following perturbed equations of motion 



5 x 


= AS x 


+ 


B6 u 


( A- 9 ) 


6 y 


= 05 x 


+ 


D6 u 


( A-10 ) 



The two equations above approximate the dynamic behavior 

of the system about the operating point (x , u ). The 

s s s s 

elements of the motion A, B, C and D are given by 

Ai j - i = 1 , 2 . . . n j = 1 , 2 . . . n 

5 xj J 



(A-7 ) 
( A- 8 ) 
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Bij 



= a xi 
5 u] 



Cij 



6 yi 
6 



Dij 



o u j 



The representation of equations (A-9) and (A-10) is given 
by the following block diagram. 




X 



ss 



Figure A-l. Block Diagram of a Linear Model. 
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